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Abstract 

Detection systems rely more and more on on-line or 
off-line comparison of detected signals with basis sig- 
nals in order to determine the characteristics of the 
impinging particles. Unfortunately, these compar- 
isons are very sensitive to the random time shifts that 
may alter the signal delivered by the detectors. We 
present two fast algebraic methods to determine the 
value of the time shift and to enhance the reliability 
of the comparison to the basis signals. 



1 Introduction 

Modern detection systems often rely on the on-line 
shape analysis of the signals delivered by the de- 
tectors. In this context, signal uncertainties due to 
noise, cross-talk or time jitter are of tremendous im- 
portance. This paper introduces new methods to ad- 
dress the latter effect. 

In grid search algorithms [1 , the detected signal is 
systematically compared to all the signals of a basis 
in order to find the best-match. The signals from the 
basis are considered as reference signals, they do not 
present any time shifts. The characteristics of the 



"event" that has generated the pulse are determined 
by searching the best-match among the basis signals. 
This technique is used in many modern applications, 
for example to determine the charge and mass of the 
detected ions [5] or the numbers of carbon and hy- 
drogen in atomic clusters [3J H]. We will illustrate 
the methods by the search of gamma-ray interaction 
locations [5] in an HPGe crystal [SI [7]. A basis of 
signals corresponding to a set of regularly spaced hit 
locations covering the crystal volume is generated us- 
ing a signal simulation code [S], or a crystal scanning 
device [S]. Locating an experimental gamma-ray in- 
teraction consists in finding the best-match between 
the delivered signal and the basis signals. This best- 
match is supposed to correspond to the point closest 
to the actual hit. When detected signals are altered 
by time shifts, this property may no longer be true. 

Usual methods to correct random time shifts con- 
sist whether in determining a given point on the sig- 
nal (e.g. the crossing of a threshold) or in comparing 
the experimental signal to versions of the reference 
signal translated by different At. In the former case, 
the comparison, relying on a single sample, is very 
sensitive to noise and would not allow determining 
time shifts smaller than the sample duration. In the 
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latter case, the computation time is cripplingly long 
for most on-line applications. 

We present in Section 2 a reliable and fast method 
to find the best-match to a signal affected by ran- 
dom time shifts among a basis of reference signals. 
This method does not give directly the value of the 
time shift. In Section 3, we present an algebraic fast 
method to both determine the time shift of the de- 
livered signal and compare it to a signal basis. 

2 Comparison of time-shifted 
signals with reference signals 

2.1 Method 

In the following, the signal not affected by a time shift 
will be noted so{t) where t is the time index (sample 
number). To reduce the size of the equations, the 
signals will also be denoted as vectors (bold lower 
cases). The signal actually delivered by the detector, 
s, is different from Sq due to random noise and time 
shift. The squared residue between the reference sig- 
nal and the corrected signal can be written as: 



(1) 



The pulse shape analysis algorithms generally con- 
sist in comparing detected signals with a set of ba- 
sis signals. The best-match criterion is usually the 
residue, but, when the data are affected by noise, 
fluctuations or uncertainties, the rigorous criterion is 
the chi-square: 



signals are translated by a random, unknown, time 
shift. Where the signal is horizontal, if the shift is 
small enough, it has no influence on the signal ampli- 
tude. On the other hand, for rapid increases or rapid 
decreases, the biases on the sample values are high. 
In fact, the shift is proportional to the signal slope 
and to the time shift: As(t) = —p{t) At as can be 
seen in Fig. [l] The denominators of the chi-square 
terms are the quadratic sums of the uncertainties due 
to the noise and to the time jitter: 



(3) 



where CAt is the standard deviation of the time shift, 
which value is supposed to be known. The estimate 
of the signal slope is given by: 



Pit) 



■so(t + l)-.so(t-l) 
2ts 



(4) 



where ts is the time width of the signal samples. 
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where cTg is the standard deviations of the uncertain- 
ties on s. The choice of the residue criterion is jus- 
tified by the fact that, most of the time, the signals 
are affected only by a white noise with a standard 
deviation identical for all samples. In this case the 
residue and the chi-square minimizations give the 
same results. However, when the uncertainties are 
not the same for all samples, the chi-square mini- 
mization must be chosen. This is the case when the 



Figure 1: The time shift At induces an amplitude 
shift As proportional to the slope p. 



2.2 Results 

The validation of the method will be performed using 
the AGATA segmented coaxial HPGe [IT signal. A 
signal basis has been generated along a cubic grid of 
2 mm step using the MGS code [8j . When a gamma- 
ray interacts with the crystal, signals are induced in 
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the segment of the detector where the hit took place 
and in the neighboring segments. To handle the infor- 
mation as a whole, we introduce the so-called meta- 
signals, which are the concatenation of all the induced 
signals. An example is given in Fig. |2] The signal 
corresponding to the hit segment is the fourth one. 
The other eight signals are those induced in the di- 
rect neighbors of the hit segment. The time shift 
being the same for all segment signals, it also ap- 
plies to meta-signals. In the following, the j^^ basis 
meta-signal will be noted s^-. The basis meta-signals 
play the role of the reference signals. The test meta- 
signals are also generated with the MGS code. They 
will be noted s. The sampling rate for basis and test 
signals corresponds to = 10 ns samples. However, 
the test signals are first generated on 1 ns bins so 
that they can be moved by time shifts smaller than 
the sample time, then grouped into 10 ns samples. 



grid search algorithm using both minimization crite- 
ria. The qualities of the results are characterized by 
the error on the location of the hit. A first test is 
performed considering a given point at the center of 
a segment, i.e. belonging to a region of the crystal 
where the sensitivity (of the signal shape to the lo- 
cation of the hit) is low. This means that the results 
obtained for this particular test signal are among the 
worst that can be expected from the crystal. The re- 
sults are shown in Fig. [3] which represent the error 
averaged over noise using the residue (red curve) and 
the chi-square (blue curve) criteria. As can be seen, 
the results for the chi-square are much better. Up 
to a time shift jitter of 5 ns no error is made on the 
location of the hit. For higher time shifts, the error 
is more than twice as large for the residue criterion. 

5| 1 
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Figure 2: Example of meta-signal obtained by the 
concatenating of the hit segment signal (fourth sig- 
nal) and its eight neighbors signals. The red signal 
corresponds to the reference, the black line represent 
the test signal (the time shift is 75 ns and the noise 
is 1% of the signal maximum). 




Figure 3: Error on the hit location as a function of 
the standard deviation of the time shift jitter. The 
signal corresponds to a point located at the center 
of a segment and is affected by a 1% noise. The 
red curve shows the average results using, as mini- 
mization criterion, the residue and the blue line the 
average results using the chi-square, Eqs. (|2]),([3|. 



To compare the results given by the two compari- 
son criteria, we have simulated signals affected by a 
1% noise and by time shifts sorted at random in a 
Gaussian distribution with standard deviation a^t- 
The location of the hit is then determined with a 



The results obtained for the whole set of signals 
are shown in Fig. [4] leading to the same conclusions. 

In spite of its simplicity, this method gives a very 
good improvement to the location of gamma-ray in- 
teractions. Moreover, the estimate of the signal stan- 
dard deviation given in Eq. ([3]), depends only on the 
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Figure 4: Error on the hit location as a function of 
the standard deviation of the time shift jitter for the 
whole set of signals with 1% noise. Upper plot: us- 
ing the residue criterion. Lower plot: using the chi- 
square criterion. 



reference basis and on the standard deviations of the 
noise and the time shift. Thus, the standard devi- 
ation distributions for each of the basis signals can 
be pre-calculated off-line. The increase in computer 
time for the on-line calculation of Eq. ^ instead of 
Eq. ([5]) is a low price to pay. 

3 Algebraic time shift correc- 
tion 

In the following, we search for the time shift At to 
impose to a reference signal Sq so that it matches the 
signal s delivered by the detector. In other words, we 
search for At such that: 

R^ = J2isoit + At)-s{t))\ (5) 
t 

is minimum. Here, we use the residue instead of the 
chi-square since, unlike in the previous Section, we 
need a criterion sensitive to the time jitter. Assuming 
that the time shift is small enough, the shifted signal 
can be approximated by: 

/ N dsn , 1 d^Sn . n 1 i9^Sn . o 

(6) 

The first three derivatives of the signal can be es- 
timated by: 

_ dso _ So(t+l)-So{t-l) 

fy-) — dt — 2i, 

_ (Pso _ So(t+2)-2so(f) + So(t-2) 

^(,\ _ a^SQ _ so(f+3)-3so(t+l)+3so(t-l)-so(t-3) 
' y^) ~ at^ ~ it J 

(7) 

The developing of the residue to the second order 
gives: 

R\At) = (so - s)2 

+ 2p.(so-s)At (8) 
+ [p2 -t- q.(so - s)] A^ . 

The minimum of this functiorQis reached for: 

^As we need /?^(At) to show a minimum, it must be devel- 
oped at least to the second order. 
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At = - 



c — p.s 



(9) 



b — q.s 

where c and b are constants that depend only on Sq 



p so , 



b = 



q.so . 



(10) 
(11) 



Thus, they can be pre-calculated. 

In order to correct longer time shifts, the third or- 
der term of the residue can also be considered. This 
term, to be added to Eq. (Isl), reads: 



-[ir.(so -s) +p.q] Ai3. 



(12) 



Once more, setting the derivative to zero gives the 
value of the time shift that minimizes the residue: 



At 



-{b - q.s) + ^J{b - q.s)2 - 2(a - r.s)(c - p.s) 

a — r.s 

(13) 

where a is also a constants that depends only on Sq: 



3p.q + r.so . 



(14) 



The quality of the match of two signals can also 
be measured using their scalar product [101 E] • The 
resulting time shift is At — — (p.s)/(q.s). The tests 
show that the agreement with the true time shift is 
much better using the minimum of the residue than 
using the maximum of the scalar product. Thus only 
the results obtained with the residue are shown in the 
next Section. 

3.1 Results 

3.1.1 Validation criteria 

In order to validate this time shift correction method, 
we consider two criteria. The first one is obviously 
the error on the time shift defined as the difference 
between the real time shift and the values found by 
Eqs. (|9| and (13 1. The second criterion measures the 
error on the location of the hit. In order to evaluate 
the improvements, we will compare the errors, with 
and without time shift correction, using grid search. 



3.1.2 Error on the time shift evaluation 

We first consider as test signal a meta-signal cor- 
responding to the center of a segment. This test 
signal was translated by 1 ns time steps from 1 to 
100 ns. Three deposited energies are considered: 300 
keV (the noise is 1% of the maximum amplitude of 
the segment signal), 100 keV (3% noise) and 30 keV 
(10% noise). 

The results obtained by the second order cor- 
rection, Eq. (|9]), for this signal are presented in 
the upper plot of Fig. [5] The time shift error 
Aitrue — Atfound IS shown as a function of the im- 
posed time shift for different amount of noise. The 
bold lines represent the mean time shift error aver- 
aged over noise and the upper and lower dashed lines 
the standard deviations. As can be seen, the results 
are very good. In this case, the time shift correction 
is better than 1 ns up to time shifts equal to 15 ns and 
the error is less than 10% up to 25 ns. Remarkably, 
the average correction is not sensitive to the amount 
of noise, even for a large value. 



The results for the third order correction, Eq. ( 13 ) 



are show in the lower plot of Fig. [s] In the left part 
of the curve, the quality of the correction is the same 
as for second order correction. As expected, the time 
shift correction is better for large values of At. 

In these plots, we have presented the results ob- 
tained in a large range of time shifts. In fact, when 
the shift is longer than some sample times, other effi- 
cient correction procedures may be used. The residue 
minimization method is particularly suited for small, 
intra-sample, adjustments. Thus, in the following, we 
will focus on the left part of the figures. 

We now consider the whole set of basis signals. 
Each of them is successively used as test signal be- 
ing added a 1% noise and shifted by 1 ns steps. The 
results of the correction are shown in Fig. [6j The 
upper plot corresponds to the second order correc- 
tion. In the first sample, that is for time shifts lower 
than 10 ns, the error on the correction is lower than 
0.25 ns in 68% of the cases. For a two-sample shift, 
the standard error of the correction is better than 
10%. Inside the first sample, the correction obtained 
with the third order method is even better than with 
the second order one. A small positive bias is visible 
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in the second sample. 
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Figure 5: Error on the time shift correction as a func- 
tion of the imposed time shift for a point located at 
the center of a segment. Upper plot: second order 
correction, Eq. Q. Lower plot: third order correc- 
tion, Eq. (13). The red line shows the result when 



no correction is applied. The bold lines correspond to 
the values averaged over noise and the dashed lines 
to one standard deviation from the mean (blue: 1% 
noise, violet: 3 % noise, green: 10% noise). 
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Figure 6: Error on the time shift correction Attruo — 
Aifound as a function of the imposed time shift for the 
whole set of signals. The black full line corresponds 
to the average value obtained for the test signals and 
the dashed lines to one standard deviation from the 
mean. The red line shows the result without correc- 
tion. Upper plot: second order correction. Lower 
plot: third order correction. 



The final goal of pulse shape analysis, in the case 
of segmented HPGc detector, is to determine the lo- 
cation of the hits. As for the chi-squarc method, we 
now test the improvement on the position. 
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3.1.3 Error on the hit location 



The correction of the time shift allows a better lo- 
cation of the gamma-ray interactions. In the case of 
Fig. [7] when no correction is applied, a 4 ns shift 
may result in a 2 mm error on the position. On the 
other hand, a second order correction allows to have 
a perfect hit location up to a time shift of 15 ns, (1 
% noise) which is largely sufficient for signal decom- 
position purpose. 



16 
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Figure 7: Error on the location of a hit as a function 
of the imposed time shift. The signal is the same as 
in Fig. |3] and |5] The red curve shows the average 
results without correction and the blue curve shows 
the average results after second order correction (the 
dotted lines correspond to the mean plus or minus 
the standard deviation). 



The results obtained for the whole set of test sig- 
nals are shown in Fig. [8] The upper plot corresponds 
to the error on the location of the hit when no time 
shift correction is applied and the lower plot present 
the result after second order time correction. When 
the correction is applied, the grid search makes al- 
most no location error when the time shifts are lower 
than the sample time. 
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Figure 8: Error on the hit location as a function of 
the imposed time shift for the whole set of signals 
with 1% noise. Upper plot: no time shift correction. 
Lower plot: second order time shift correction. 
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4 Conclusions 



References 



This paper proposes two different methods for the 
correction of signal time shifts. The robustness of 
these methods is connected to the fact that, unlike 
CFD threshold techniques for example, they use the 
whole signal to evaluate the shift, thus the amount 
of relevant information is high and the negative in- 
fluence of noise is minimized. 

The first method allows the comparison of basis 
signals with detected signals affected by white noise 
and random time shift. It improves the confidence 
on the determination of the best-match basis signal. 
In our example, this induces a much better location 
of the gamma-ray interactions in a germanium de- 
tector. The increase of computing time, due to the 
use of a chi-square rather than a residue, is largely 
smaller than what would have been necessary to make 
a minimization by systematically comparing the de- 
tected signal with a set of time shifted versions of the 
reference signals. 

In the second part of this paper, we have proposed 
an algebraic evaluation of the time shift based on the 
Taylor polynomial of the reference signal and on the 
minimization of the residue. The reliability of this 
method is comparable with other methods when the 
shift is greater than a few samples. For intra-sample 
time shift, the results are excellent (the standard de- 
viation is lower than one fortieth of the sample time 
for a 1% noise) . The choice between second and third 
order corrections depends on the application and on 
the available computing time in case of on-line cor- 
rection. The second order correction entails mostly 
the calculation of two scalar products, Eq. ([9|, and 
the third order correction, the calculation of three 



scalar products and a square root, Eq. (13). In prin 



ciple, third order correction gives better time shift 
estimates, but, as it entails the calculation of higher 
degrees functions, the uncertainty on the result is 
larger. However, whatever the order of the correc- 
tion, the method is much faster than the traditional 
method of exhaustive time-shift trials. 
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